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We consider the non-minimally kinetic coupled version of DGP brane model, where the kinetic 
term of the scalar held is coupled to the metric and Einstein tensor on the brane by a coupling 
constant £. We obtain the corresponding held equations, using the Friedmann-Robertson-Walker 
metric and the perfect huid, and study the inhationary scenario to confront the numerical analysis 
of six typical scalar held potentials with the current observational results. We hnd that among 
the suggested potentials and coupling constants, subject to the e-folding N = 60, the potentials 
V((j>) = cr0, V(4>) — erf 2 and V(f) = erf 3 provide the best hts with both Planck+WP+highL data 
and Planck+WP+highL+BICEP2 data. 
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I. INTRODUCTION 

During the last several years, the braneworld scenario has been considerably studied within the variety of different 
models. According to this scenario, we are living over a three-dimensional hypersurface in a higher-dimensional 
spacetime; the standard model particles are conhned on the brane, and the gravitons propagate in the bulk spacetime 
m- Using this scenario in a hve-dimensional spacetime, Randall and Sundrum (RS) proposed two kinds of resolutions 
for the hierarchy problem [3]. They showed that in a five-dimensional spacetime one may derive the effective Einstein 
equations for the 4 D brane metric obtained by projecting the 5 D metric onto the braneworld which result in the 
most general form of the 4 D gravitational field equations for a braneworld observer gnu. Induced gravity brane 
model proposed by Dvali, Gabadadze and Porrati (DGP) is another example of this scenario to account for the self 
accelerating behaviour of the universe [8j. Many authors have studied the geometrical mm as well as the cosmological 
mmm aspects of this new 5 D gravitational model. DGP model with a bulk cosmological constant and a tension 
of the brane, with energy scale much larger than the 5 D Planck mass, leads to the effective cosmological constant on 
the brane which is extremely reduced in contrast to the RS model, even if the cosmological constant and the tension 
are not fine-tuned [18] . 

The inflationary scenario can resolve the problems of standard cosmology such as the flatness, horizon, monopole 
and relics problems. In most of the successful inflationary models the universe is filled with a scalar field so called 
Inflaton whose potential energy is dominant over its kinetic energy [T9H26] . However, several problems remain without 
concrete solutions |24kF26l 128]. Hence, many other inflationary models such as the braneworld models [291133] . models 
with non-minimally coupled inflaton field [f34H4T| . modified gravity [42, [44] [45] . and models with a wide range of 
potentials have attracted so much attention in the recent years. In this regard, variety of models have been proposed 
however those models are viable that show consistency with observational data and provide us with a mechanism 
for generating the initial fluctuations and perturbations in the early universe as the seeds for the formation of the 
structures in the universe. In such models, the fluctuations in the scalar field as well as the transverse and traceless 
parts of the metric lead to the scalar and tensor power spectrum, respectively [T91I26] . The scalar power spectrum 
is nearly scale-invariant, with the order of unity, and the good point is that the exact value of spectral index can be 
obtained by using the observational data. Moreover, the running of spectral index and the tensor-to-scalar ratio can 
also be constrained observationally. Comparison between the calculated values of these parameters and the recent 
observational data are the most powerful probes for ruling out or keeping a specific inflation model. 

Such a study has already been done in the context of non-minimal DGP braneworld inflation in Ref. m, where 
the non-minimal feature of the model was attributed to the non-minimal coupling between the inflaton field and the 
induced Ricci scalar on the brane. Also, the observational constraint were analyzed with respect to the background 
of Planck+WMAP9+BAO data and the potential V(f) ~ exp (—/3f) was obtained as the best fit case. 

Here, we develop a study similar to Ref. m to find other possible inflation models, consistent with the observations, 
in DGP braneworld scenario. However, the present study is much different from Ref. (27' in four senses. The first is 
that here we follow a different approach to use the non-minimal feature in our model. Rather than considering the 
non-minimal coupling between the inflaton field and the induced Ricci scalar on the brane, we consider a non-minimal 
coupling between the kinetic term of the inflaton field and Einstein tensor on the brane. Such models are known as 
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non-minimal kinetic coupled gravity [48] . The second is that here we analyze our observational constraint with respect 
to the background data of Planck+WP+highL+BICEP2 rather than Planck+WMAP9+BAO. The third is that here 
we consider six types of inflaton potentials, more or less different from those of Ref. [27]. and perform a numerical 
analysis on the inflationary parameters of this model to confront them with Planck+WP+highL+BICEP2 data. The 
forth is that here we obtain three new best fit potentials rather than (and different from) one obtained in Ref. [27]. 
It was already found that some potentials which are suitable for inflation in 4-dimensional model, cannot lead to 
a successful inflation in the minimal case of 5-dimensional DGP model. Moreover, some potentials which are not 
compatible with observational data in a 4-dimensional model, can lead to viable results in a minimal 5-dimensional 
DGP model [27]. In this paper, by considering a non-minimal coupling between the kinetic term of the scalar field 
and Einstein tensor in a 5-dimensional DGP model, we obtain new scalar field potentials suitable for inflation. 


II. FIELD EQUATIONS IN THE BRANE SCENARIO 


We assume a 5 D bulk spacetime (A4, ^ Qab ) with the coordinates X A (A = 0,1, 2, 3, 5) and a 4 D brane (M,g^ u ) 
located at a hypersurface B(X A ) = 0. The standard action for the braneworld is written as 


S = ^Bulk + 5br 


(1) 


Sbu lk= [ d 5 Xv /Z ^[A (5) i?+ (5) L: 

JM 


( 2 ) 


*^brane 


[ d A xJ^g 

JM 


2^ T T^rane (g a /3 , 'Ip) 


(3) 


where r\ corresponds to the 5 D gravitational constant ms, and ^R are the 5 D the matter Lagrangian in the 

bulk and scalar curvature, respectively. Also, x v (v = 0,1,2,3) are the induced 4 D coordinates on the brane, K± is 
the trace of extrinsic curvature on either side of the brane IEBEU and Tbran e(gaPi'ip) is the effective 4 D Lagrangian 
which is given by a typical functional of the brane metric g a p and matter fields ip . 

The five-dimensional Einstein equations in the bulk are given by 


^Gab = ^ ^Tab + t ab 8(B) 


where 


G) 


Tab = -2- 




S(5)gA B 

is the energy-momentum tensor of bulk matter fields, and 


(5) 


(5) L n 


9AB , 


^-^brane , ^ T 
^~jjV = - ~ — T C/^vL /brane : 


Sg^ 


(4) 


(5) 


(6) 


is the effective energy-momentum tensor localized by 5(B) on the brane. The induced 4 D metric can be written as 
g_\B = ^9ab ~ nA’iiB, where n .\ is the spacelike unit-vector field normal to the brane hypersurface M. Following [3], 
0, and m one obtains the gravitational field equations on the braneworld as [? ] 


r _ H 

,1V ~ 3 


{5) T rs 9*9 S „+9^( (5) T rs n R n s - 1 (5) T 


D v t" ^-2^T RS n R g\, 


^5 ixv , 


(7) 

(8) 


where 


and 


a/3 


^~fiaTif T -j 2 TT I 1U 


~rp {h) r M N R S 

= j Cmrns n n 9 u 9 


24 


9^r 


H U v' 


(9) 


( 10 ) 













3 


III. DGP BRANE’S MODEL WITH NON-MINIMAL KINETIC COUPLED GRAVITY 

In this section, we modify DGP braneworld model by a non-minimal kinetic coupling term in the Lagrangian 

ibrane = ^ R ~\ W" ~ C G"") “ V(<f>) ~ A + L m , (11) 

where /i is a mass scale which may correspond to the four dimensional Planck mass 777 , 4 , R is the Ricci scalar, ( is a 
coupling parameter with dimension of ( length ) 2 , V{(j)) is the scalar field potential, A is the tension of the brane, and 
L m is the Lagrangian of other matters on the brane. The presence of Einstein tensor in the kinetic term of the inflaton 
field is novel and casts this model in the context of non-minimal coupled gravity. Also, we take only a cosmological 
constant in the bulk. 


A. Field Equations on the brane 

To obtain the field equations on the brane, we calculate the energy-momentum tensor of the brane as 

< = -A^ + n - ^ + W + C©£, 

where 


W = - -S^v p (t>v p <f) - KV{4>), 


O/zzy — — + 2V a 0 

4(V0) 2 G^ + ^[ - § V a V a V^ + \{U^) 2 - V a 0 V^iT^]. 


Substituting this equation into Eq.([7|), 
A 


where 


one can find the effective equations for 4 D metric as [18] 

A 


1 + — + /^5/i 2 /C M ^ pcr (T a ^)G pcr + A 




■E, 


'pu-> 


fcpupa — ^ (dpvTpcr QppTycr Qua Tpp) + 12 + T (g^ p g ua QpuQpa) 




( 5 )a + -4a 2 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


fP = -\T m T« + + \ 9llv T a ^ - ± 9liV T 2 , 


pu 


(19) 


4P — -^Gp a G u a + —GGp V + -gp, u G a pG al3 - —g^G 2 , 


( 20 ) 


T and G being the trace of energy-momentum and Einstein tensors, respectively. Because of the Bianchi identity, the 
Codazzi equation reads as D v Tp V = 0 which implies the energy momentum conservation, i.e. 


D u Tp V = 0 . 


( 21 ) 
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B. Cosmology of non-minimal kinetic coupled DGP model 


We take the spatially flat isotropic and homogeneous FRW line element on the brane 

ds 2 = —dt 2 + a 2 (t)dijdx l dx^ , 


( 22 ) 


where Sij is a symmetric 3-dimensional metric and a(t) is the scale factor. Studying such universe with a perfect fluid 
and following [4], we can write = 0 implying that 


= 0 . 


as 


The original field equations 0 can be written 

g°„ = "„-n. 


where 


and 


with 


G'i = 


G° 0 = -3 ( H 2 + ) , 


G'j = -(2E + 3fC + ^)S\,, 




71-4 j = ( T °0 - 2r \) 


r° 0 — —(A + p) — ^ 2 G° 0 + -(1 + 9 (H 2 )<fi 2 — V(4>), 


(23) 

(24) 

(25) 


(26) 


(27) 


(28) 


t \ = (P~ A )S i j - n 2 GP + i2 


C(2 


1 + C ( 2H + 3H Z + 


4 W -1 )] S*j - V(</>)6\ 


Eqs.(24) and (25) are written as 


3X = i( 5 )A + E° 0 + (-A - p + 3v 2 X + ^(1 + 9C77 2 )0 2 - V(0)) 


(29) 


(30) 


1 + ^-A-p + 3/l+- ( l + PV-^) (^ 2 + | 


42 


y 


= -| E °o + if (p + P + ^ 2 - 2 QHU + j^ 2 T) (—A - p + 2>n 2 X + 1(1 + gC^ 2 )^ 2 - F(0)) , (31) 


where 


X = H‘* + 

a z 

Y = H- 4. 


(32) 


By using Eq. (23), we can obtain the equation of motion for £’° 0 as follows 

£° 0 + 4HE° 0 = 0. 


(33) 
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By integrating from the above equation we can easily find 


E 


£o 


( 34 ) 


where So is an integration constant. Now, we must solve the equation (30), as a quadratic equation with respect to 
X , which can be written as 



1 

v 


Pm + Po (l + eA(p, a)) 


where 


Pm — P + P(f> — p + - (1 + 9 (H 2 )(f) 2 + V ((/)), 


and 


Po = m\ + 6^- , 

/i 2 

with the mass scale tua = A 1 / 4 . Also, e stands for either +1 or —1, and *4 is defined by CHI 


(35) 


(36) 


(37) 


where 


A = 



{^Prri p 



1 

2 


A 0 = — , 

V Po 

r = 6 (o < r < i). 

PoP 2 


(38) 


(39) 

(40) 


Equation ( [35] ) is considered as the Friedmann equation of this model. Note that the choice of sign for e has a 
geometrical meaning m and it is determined by the initial condition of the universe. 


IV. INFLATION 


In the slow-roll regime, we have 9 (H 2 <j) 2 V ((/)), (j) 3 H<fi [15| and because of (34) at inflationary stage (a So), 

we may ignore the integration constant So by setting So = 0. So, the energy density takes the following form 


V (</>) = V. 


(41) 


By varying the Lagrangian © with respect to the scalar field we have 

(<j> + 3 H4>) + K{H 2 4> + 2 HH<fr + 3 H 3 (j)) + V' = 0. 

Using the slow-roll approximation, we obtain the following equation of motion 

3H<j> (l + C(2 H + 3if 2 )) + V' = 0. 

So, the Einstein equations in slow-roll approximation take the following forms, respectively as 


H 2 = 


3 p 2 


V (<f>) + po (l + eA) 


H = - 


V 


a 


V 9H 2 (l + ({2H + 3H 2 )) 


1 -I - eTA 


(42) 


(43) 


(44) 


(45) 
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where 


A 


.9 V 

A 2 q + 2T — 
Po 


(46) 


and ' denotes The slow-roll parameters defined by £ = — jpi and rj = £ — take the following forms, respectively 
as 


6 = 


/i 2 U' 2 1 

[l + C(2i7 + 3 H 2 )\ 


1 + eTA 
1# ^(l + eA) 


(47) 


and 


6 — 


T] = 


6tf 2 (i+c(27r+3tf 2 )) 




(. 


3C sH 2 


\ 


\^1+C(2ij+3H 2 ) J 


1 - 


( 


\ 


yi+C(2H+3H2) J 


(48) 


where 

• 4 = (• 4 S + 2 r D '' 

B = l + eIVt _ \ 

C = l + y(l + eA). (49) 


The number of e-folding is given by N = fj f Hdt where (t t ) and (t e ) are the initial and end time of inflation, 
respectively. For a warped DGP model with a non-minimally kinetic coupled gravity on the brane, we will get the 
following expression 


N = — 


1 

/i 2 




1 + C(2 H + 3 H 2 ) 



(50) 


where pi and <pf are the values of <p when the radius of universe crosses the Hubble horizon during inflation and 
exits the inflationary phase, respectively. A useful tool to test the viability of inflationary models is the spectrum 
of perturbations produced due to the quantum fluctuations around their homogeneous background values. The 
conformal-Newtonian metric is given by [571459] 


ds 2 = a 2 (rj) 


— (l + 2^/) dp 2 + (l — 2^)Sij dx l dx J 


where T is called the Bardeen potential. The primordial power spectrum is defined by the following expression 

t2 \ 2 


P R = (K Y 

V2t tc!>) 


(51) 


(52) 


k~aH 


Using the scalar field equation of motion in slow-roll regime (i.e. equation(|43|)), we get 


p„. (*\ 

V2ir 4.) 


V' V 


k=aH 


12tt V V V ' 


i 2 


1 + C(2 H + 3 H 2 ) 


1 + ^(1 + ^) 


1 3 


(53) 


In the slow-roll regime, we know that e <C 1 and H ~ 0 ; so the equation (53) takes the following approximate form 

/ \ 3 

V 


P„ -- 


V 


12tt 2 /x 6 l V' 


(l + 3C^ 2 ) (l + (l + eA) ) . 


(54) 
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The scalar spectral index, which describes the scale-dependence of the perturbations is defined as 

dlnP 2 

n a — 1 = 


dink 


(55) 


where dlnk(q !>) = dN((/>). As is seen, for n s = 1 the power spectrum of the perturbation is scale invariant. In our 
warped DGP model, we obtain the scalar spectral index in the slow-roll regime as follows 


_ /r/^v 

1 “ - ( V ) 


2 f 3(1 + el\4 -1 ) 2V"V 


l + $(l + eA) 


V ' 2 


(\ + Q{2H + ZH 2 ))(l+ P ^(\ + eA)) 


(56) 


where A is defined in (49). The running of spectral index in our model is given by 
dlnn o 


a = 


dink 


— 6 


+ 2 


[l + C,(2H + 3H 2 ]C 

v \ 2 , y 




' 2 ,V\ 2 , 2 V \ Jn 2 (V'\ 

" (y> x bV )- 4 Mv) 


1 + C(2 H 

/ o\ 2-| 


Y 


c 

TT + 


CCF 


3H 2 )\C ' M 2(i 

b 2 £) 




C(2H + 3H 2 )) y V C 

(1 + C(2^ + 3^2))C \2\VJ ) V C 2 v 1 



(57) 


where 


.4 


= Ul 


u V 
2T — 
po/ 


= l 


eUM, 

C = 1 + (l + e^4) . (58) 

The tensor perturbation (gravitational wave) amplitude of a given mode, at the time of Hubble crossing, is given by 

2 




/i 2 \2tt 


k=aH 


In our model and in the slow-roll regime, we find the primordial tensor perturbation 

8 V 


Pi = 


(27r) 2 3/i 2 


( 1 + 7 ( 1 +^)). 


and the tensor spectral index is given by 


dlnPl 


n„ = 


dink 

so, in the slow-roll regime we can express it as follows 

n T = —2s . 

Now, we evaluate the tensor-to-scalar ratio as 




PI 


fP 

2 


V'\ 2 


(v) 


16 


(1 + ((2H + 3H 2 )) 11 + f (1 + eA) ] 2 


(59) 


(60) 


(61) 


(62) 


(63) 


Up to now, we have presented the equations of cosmological dynamics. In the following, we perform numerical 
analysis on the inflationary parameters of the warped DGP model with a non-minimally kinetic coupled gravity on 
the brane. Now, we shall consider some types of potentials by substituting them in the integral of equation (50) and 
then solve this equation. But, first we should find the value of </> at the end of inflation, namely 0/, by setting s = 1 
in Eq.(|47|) (corresponding to the end of inflation). Then, we put it in (50) and find fa in terms of N and substitute 
fa in n s , r and cq for any given values of N. Now, we can compare these inflationary parameters with the recent 
observational data. 
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FIG. 1: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ = 10 -6 and the potential V(<f>) = The plots are depicted for six 

values of number of e-folding N. For all given viable values of N , the non-minimally kinetic coupled DGP model lies inside 
the 95% CL of the Planck+WP+highL data, but it does not lie in the 95% CL of the Planck+WP+highL+BICEP2 data. 
The values of running of spectral index are very close to zero for all given values of N , and coincides with the prediction of 
the single monomial chaotic inflationary models with 50 < N < 60 (For comparison, there is a purple strip hidden behind the 
running spectral indexes of non-minimally kinetic coupled DGP model, showing the running spectral index of single monomial 
chaotic inflationary models with 50 < N < 60). The non-minimally kinetic coupled DGP model lies inside the 95% CL of the 
Planck+WP+BAO:ACDM dn s /dlrik data, but it does not lie in the 95% CL of the Planck+WP+BAO:ACDM dn s /dlri } c + r 
data. 


V. OBSERVATIONAL CONSTRAINT 

In this section, first we introduce our model parameters as 

~ ra^ -3 , fi ~ 777/4 t A ~ 0 , < 5 >A ~ 0 , T - 0.99 , k ~ 0.002Mpc -1 , (64) 

together with two parameters ((, a); the first one is the nonminimal coupling constant for which we shall take different 
values so that one can decide which one shows best fit with observations and the second one is the parameter defined 
in the scalar field potential (se bellow) which, with no loss of generality, we will take its value to be of the order of 
unity. Then, we investigate these models and compare the results with the observational data. 


A. V ((/)) = a<j) 2 

The non-minimally kinetic coupled DGP model is well inside the joint 95% CL Planck+WP+highL data (red area) 
for all values of e-folding TV, but it does not lie inside the joint 95% CL Planck+WP+highL+BICEP2 data (blue 
area). In the left plot of Fig|l] the behavior of tensor to scalar ratio versus scalar spectral index is shown in the 
background of the Planck+WP+highL+BICEP2 data for six values of N . In the right plot, the evolution of running 
of spectral index versus scalar spectral index has been plotted for the similar situation. It is seen that, for all six 
values of the number of e-folding, the running of scalar spectral index is negative and close to zero. 


B. V{<t>) = a<j>i 

The non-minimally kinetic coupled DGP model is well inside the 95% CL of the Planck+WP+highL data, but it 
does not lie inside the 95% CL of the Planck+WP+highL+BICEP2 data. Evolution of tensor to scalar ratio versus 
scalar spectral index is shown in the left plot of Figj2j For all given values of V, the non-minimally kinetic coupled 
DGP model lies inside the 95% CL Planck+WP+highL data. Evolution of the running of scalar spectral index versus 
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FIG. 2: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a non- 
minimally kinetic coupled DGP model with £ = 10 -6 and V(4>) = cr0 2//3 potential. The figure has been plotted for six values of 
the number of e-folding. For all given viable values of TV, the non-minimally kinetic coupled DGP model is well inside the 95% 
CL Planck+WP+highL data, but it does not lie within the 95% CL Planck+WP+highL+BICEP2 data. The values of running 
of spectral index are very close to zero for all given values of N. The non-minimally kinetic coupled DGP model lies inside the 
95% CL of the Planck+WP+BAO:ACDM dn s /dlrik data, but it does not lie in the 95% CL of the Planck+WP+BAO:ACDM 
dris/dlrik + r data. 


scalar spectral index has been plotted in the right panel of Figj2] For this case, the running is negative and close to 
zero. 


C. V(0) = (7 (f> 

A minimally coupled four-dimensional setup with this potential lies within the 95% CL of the 
Planck-PWMAP9+BAO data [26]. Our braneworld model (non-minimally kinetic coupled model), with this lin¬ 
ear potential, lies within the 95% CL Planck+WP+highL+BICEP2. A minimally coupled DGP model with this 
potential lies still inside the 95% CL of the Planck+WMAP9+BAO data. As before, we consider six values of number 
of e-folding. In the left plot of Fig|3j we see the evolution of tensor to scalar ratio versus scalar spectral index. From 
our numerical analysis it appears that in a DGP model with non-minimally kinetic coupled gravity, the model lies in 
the 95% CL Planck+WP+highL+BICEP2, for all given values of N. The right plot of Fig^] shows the evolution of 
running of scalar spectral index versus scalar spectral index. As the figure shows, for a non-minimally kinetic coupled 
DGP model with a linear potential, the running of scalar spectral index is close to zero. 


d. y(0) = cr0 2 

In EZ], it has been shown that in 4-dimensions the model with this potential lies outside and inside the 95% CL 
of the joint Planck+WMAP9+BAO data for N = 50 and N = 60, respectively. Now, we explore the situation for a 
5-dimensional model. According to the WMAP7+BAO+Ho data m, a warped DGP model with minimally coupled 
scalar field and with a squared potential, lies inside the 95% CL for TV < 70. Now, with recent BICEP2 date, the 
situations change considerably. In a minimally coupled DGP model with a quadratic potential, for all N > 40, the 
model is outside the joint 95% CL of the Planck+WMAP9+BAO data. In our model for a non-minimally kinetic cou¬ 
pled DGP model, for all given values of N the model is well inside the joint 95% CL Planck+WP+highL+BICEP2 
data. The left plot of Figj4] shows the behavior of tensor to scalar ratio versus scalar spectral index in the back¬ 
ground of the Planck+WP+highL+BICEP2 data. This figure has been plotted for six values of N. Also, we 
have plotted the evolution of running of scalar spectral index versus scalar spectral index in the background of 
the Planck+WP+highL+BICEP2 data (the left panel of Fig.4). We see that, for all six values of the number of 
e-folding, the running of scalar spectral index is close to zero. 
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FIG. 3: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ = 10 -6 and a linear potential V(<f>) ~ 0. The figure has been plotted for six 
values of number of e-folding. For all given viable values of N , the non-minimally kinetic coupled DGP model is well inside the 
95% CL Planck+WP+highL data and is almost inside the 95% CL Planck+WP+highL+BICEP2 data. The values of running 
of scalar spectral index are almost zero for all given values of N. The non-minimally kinetic coupled DGP model lies inside the 
95% CL of the Planck+WP+BAO:ACDM dn s /dlrik data, but it does not lie in the 95% CL of the Planck+WP+BAO:ACDM 
dn s /dlrik + r data. 
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FIG. 4: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ = 10 -6 and a quadratic potential V((j>) (j ) 2 . The figure has been plotted 

for six values of number of e-folding. For all given viable values of N , the non-minimally kinetic coupled DGP model is well 
inside the 95% CL Planck+WP+highL+BICEP2 data and is almost inside the 95% CL Planck+WP+highL data. The values 
of the running of spectral index are almost close to zero for all given values of N. The non-minimally kinetic coupled DGP 
model lies almost inside the 95% CL of the Planck+WP+BAO:ACDM dn s /dlrik data, but it does not lie in the 95% CL of the 
Planck+WP+BAO:ACDM dn s /dlrik + r data. 


E. V(</>) = c 74 > 3 

It has been shown in m that in natural inflation this potential lies in the 95% CL Planck+WP+highL+BICEP2 
and also it has been confirmed with WMAP9 m and Planck m data that a model with a cube potentia in 4- 
dimensions lies outside the 95% CL. In our branworld model, we obtain a different result: a non-minimally kinetic 
coupled DGP model with this potential, lies inside the 95% CL Planck+WP+highL+BICEP2 data for each given 
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FIG. 5: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ = 10 -6 and a cube potential V(<j>) ~ £> 3 . The figure has been plotted for 
N = 40, 50, 60, 70, 80 and 90. For all given viable values of AT the non-minimally kinetic coupled DGP model is well inside 
the 95% CL Planck+WP+highL+BICEP2 data and is almost inside the 95% CL Planck+WP+highL data. The values of 
running of spectral index are approximately zero for all given values of AT. The non-minimally kinetic coupled DGP model lies 
approximately close to the 95% CL of the Planck+WP+BAO:ACDM dn s /dlrik data, but it does not lie in the 95% CL of the 
Planck+WP+BAO:ACDM dn s /dlrik + r data. 


value of AT. The results are shown in Figj5] Note that the evolution of running of scalar spectral index corresponding 
to the cube potential is shown in the right plot of Fig|5] The value of running of scalar spectral index is close to zero. 


F. V(</>) = cr cj) 4 

A minimally coupled 4-dimensional model with this potential lies within the 95% CL of the Planck+WMAP9+BAO 
data m- For AT < 70, our braneworld model (non-minimal kinetic coupled) with this potential lies inside the 95% CL 
Planck+WP+highL+BICEP2 data. But for all given values of AT the model lies inside the 95% CL Planck+WP+highL 
data. The running of spectral index is plotted in the right plot of Fig.6. In the left plot of Fig.6, we see the evolution 
of tensor to scalar ratio versus scalar spectral index. 


G. V(</>) = a4> 5 

Similar to the other cases, we consider six values of the number of e-folding. In the left plot of Fig.7, we see the 
evolution of tensor to scalar ratio versus scalar spectral index. For all given values of AT, the non-minimally kinetic 
coupled DGP braneworld model with this potential lies inside the 95% CL Planck+WP+highL data, but does not lie 
inside the 95% CL Planck+WP+highL+BICEP2 . The running of spectral index is plotted in the right plot of Fig.7 
and it is close to zero. 

Our numerical analysis for a DGP model with a non-minimally kinetic coupled provides us with the value of each of 
the parameters & (the value of (f) at the beginning of inflation), </> e (the value of (f) at the end of inflation), Hi (the value 
of H at the beginning of inflation), H e (the value of H at the end of inflation), £,77 (the slow-roll parameters), V((/>i) 

(the value of potential at the beginning of inflation), and 9(to control the slow-roll condition 9£JT 2 4 > 2 « vm 
for AT = 60 (the number of e-folding), in the following tables. These values can tell us “which potentials are in best 
agreement with the recent observations.” 

It turns out that the potentials V((f)) oc </> 2 and V((f)) oc (f> 3 provide respectively the best fits with the observations 
(see Fig.4, Fig.5 and table.II). 
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FIG. 6: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ = 10 -6 and a potential of type V(<j>) oc 0 4 . The figure has been plotted for 
N = 40, 50, 60 70, 80, and 90. For all given values of N the model lies inside the 95% CL Planck+WP+highL data and lies 
inside the 95% CL Planck+WP+highL+BICEP2 data just for N = 40, 50, 60, 70. The running of spectral index is plotted in 
the right panel and it is close to zero. Neither of N = 40, 50, 60, 70 lies inside the 95% CL of the Planck+WP+BAO:ACDM 
dn s /dlrik data and the Planck+WP+BAO:ACDM dn s /dlrik + r data. 
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FIG. 7: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ =t 10 -6 and a potential of the type V(<j>) ~ 0 5 . The figure has been plotted 
for N = 40, 50, 60, 70, 80, and 90. For all given viable values of iV, the non-minimally kinetic coupled DGP model lies inside 
the 95% CL Planck+WP+highL data, but it does not lie inside the 95% CL Planck+WP+highL+BICEP2 data. The running 
of spectral index is plotted in the right plot of the figure; it is negative for iV > 50, and is positive (close to zero) for N = 40. 
The model with values of N = 70, 80, 90 and N = 80, 90 lie inside the 95% CL of the Planck+WP+BAO:ACDM dn s /dln \ c 
data and the Planck+WP+BAO:ACDM dn s /dlrik + r data, respectively. 


VI. CONCLUSION AND REMARKS 

In this paper, we have considered a 5 D bulk spacetime together with a single 4 D brane and derived the effective 
4 D gravitational equations. Then, we have studied the non-minimally kinetic coupled version of a braneworld gravity 
proposed by Dvali, Gabadadze, and Porrati, so called DGP model. We have derived the field equations, using the 
FRW metric accompanied by the perfect fluid, and studied the inflationary scenario in this model. Finally, we have 
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TABLE I: 



<J 

c 

N 

4>e 

H e 

<t>i 

H i 

£ 

V 

~9 JhJJP — 

V(0i) 

1 

crcj) 2 

1 

0.001 

60 

0.3535 

0.25705006\/3 

7.76 

0.5564\/3 

0.0020 

-0.0062 

0.007755 

2.7845 

2 

crcj) 3 

1 

0.001 

60 

0.4113 

0.2 5943 584 \/3 

8.942 

0.6919\/3 

0.0027 

-0.0054 

0.01856 

4.308 

a<}> 

1 

0.001 

60 

0.7069 

0.2 80 30475 \/3 

10.977 

1.1024\/3 

0.0041 

-0.0040 

0.11967 

10.939 

crcjp 

1 

0.001 

60 

1.4136 

0.47123562 \/3 

14.782 

4.9273\/3 

0.0075 

0.000001 

47.7464 

218.5075 

crcj)^ 

1 

0.001 

60 

0.00008 

0.000 0480 3 \/3 

13.514 

16.5601\/3 

0.0071 

0.00098 

6091.7848 

2468.1470 

c t4 d 

1 

0.001 

60 

0.0017 

0.00004805^ 

11.047 

40.6787-^/3 

0.0041 

0.0007 

2.218 X 10 5 

14892.8363 

crcj )^ 

1 

0.001 

60 

0.0080 

0.00004811^ 

7.902 

58.5089-^/3 

0.0063 

0.0012 

9.49 X 10 6 

30809.5338 

crcj)^ 

1 

0.001 

60 

0.0203 

0.00004820^ 

6.731 

101.6523^ 

0.0042 

0.0094 

8.46 X 10 6 

92998.8 

aet 

1 

0.001 

- 

” Imaginary” 

- 

- 

- 

- 

- 

- 

- 



FIG. 8: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ = 10 -3 and a potential of the type V(<j>) ~ (j) n . The figure has been plotted 
for N = 60. For all given values of n, the non-minimally kinetic coupled DGP model lies inside the 95% CL Planck+WP+highL 
data, but it does not lie inside the 95% CL Planck+WP+highL+BICEP2 data, except n— 1,2. The running of spectral index 
is plotted in the right plot of the figure. Neither of the values of n lies inside the 95% CL of the Planck+WP+BAO:ACDM 
dris/dlrik + r data and all the values of n except n — 5,6 lie inside the 95% CL of the Planck+WP+BAO:ACDM dn s /dlrik 
data. 


TABLE II: 
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FIG. 9: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for a 
non-minimally kinetic coupled DGP model with £ = 10 -8 and a potential of the type V(<f>) ~ (j) n . The figure has been plotted 
for N = 60. For all given values of n, the non-minimally kinetic coupled DGP model lies inside the 95% CL Planck+WP+highL 
data except n = 3,4, and lies inside the 95% CL Planck+WP+highL+BICEP2 data, except for n = 1/2, 2/3,6. The run¬ 
ning of spectral index is plotted in the right plot of the figure. Neither of the values of n lies inside the 95% CL of the 
Planck+WP+BAO:ACDM dn s /dlrik + r data except n — 6, and just the values n = 1,2/3,1/2 lie inside the 95% CL of the 
Planck+WP+BAO:ACDM dn s /dlrik data. 
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FIG. 10: Plots of tensor to scalar ratio (left) and running of scalar spectral index (right) versus scalar spectral index, for 
a non-minimally kinetic coupled DGP model with ( = 10 -10 and a potential of the type V{(j>) ~ (j) n . The figure has 
been plotted for N = 60. For all given values of n, the non-minimally kinetic coupled DGP model lies inside the 95% 
CL Planck+WP+highL data except n — 3,4, 5,6, and lies inside the 95% CL Planck+WP+highL+BICEP2 data, except 
n — 1/2, 2/3. The running of spectral index is plotted in the right plot of the figure. Neither of the values of n lies inside the 
95% CL of the Planck+WP+BAO:ACDM dn s /dlrik + r data, and just the values n — 1, 2/3,1/2 lie inside the 95% CL of the 
Planck+WP+BAO:ACDM dn s /dlrik data. 


confronted the numerical analysis of six typical scalar field potentials with the observational data, and found that: 

• For V(<p) = (T(j)^ and V((j>) = crcfri and the given values of TV, the non-minimally kinetic coupled DGP 
model is well inside the 95% CL of the Planck+WP+highL data, but does not lie in the 95% CL of the 
Planck+WP+highL+BICEP2. So, these potentials cannot provide the best fits with the current observations 
(see Fig.l and Fig.2). 
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TABLE III: 
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TABLE IV: 
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• For V{(j)) = (7cf) and the given values of TV, the non-minimally kinetic coupled DGP model is well inside the 95% 
CL Planck+WP+highL-f BICEP2 data. But, the evolution of tensor to scalar ratio versus scalar spectral index 
cannot provide the best fits with the current observations (see Fig.3). 

• For V{(j>) = <J(j) 2 and the given values of AT, the non-minimally kinetic coupled DGP model is well inside the 
95% CL Planck+WP+highL+BICEP2 data and the evolution of tensor to scalar ratio versus scalar spectral 
index provides the best fits with the current observations (see Fig.4). 

• For V{(j>) = ct 0 3 and the given values of AT, the non-minimally kinetic coupled DGP model is well inside the 
95% CL Planck+WP+highL+BICEP2 data and the evolution of tensor to scalar ratio versus scalar spectral 
index provides the best fits with the current observations (see Fig.5). 

• For V((/)) = cr0 4 and AT < 60, the non-minimally kinetic coupled DGP model is well inside the 95% CL 
Planck+WP+highL+BICEP2 data. Since the number of e-folding should be usually lager than 60 and because 
this potential cannot satisfy the slow-roll condition (i.e. 9 (H 2 ^ 2 <C V((/>)), it is not a good potential for inflation 
in this model (see Table II). 

• For V((p) = (J0 5 and the potentials with powers more than 5, one can show that for the given values of AT 
the non-minimally kinetic coupled DGP model is well inside the 95% CL of the Planck+WP+highL data, but 
does not he in the 95% CL of the Planck+WP+highL+BICEP2. Moreover, these potentials cannot satisfy the 
slow-roll condition (i.e 9 (H 2 <j) 2 <C V((/>)), hence cannot be considered as good potentials for inflation in this 
model (see Fig7 and Table II). 

• For V (0) = a In(0) and V (0) = ae ±( ^ and the given values of AT, we get the imaginary value of <fi at the end of 
inflation (i.e. So, these potentials cannot be considered as good potentials for inflation in this model (see 
Table II). 

• For given scalar field potentials with n = 1/2, 2/3,1, 2, 3,4, 5, 6, AT = 60 and £ = 10 -3 the non-minimally 
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kinetic coupled DGP model lies inside the 95% CL Planck+WP+highL data, and lies inside the 95% CL 
Planck+WP+highL+BICEP2 data, for n = 1,2 (see Fig.8 and table I). 

• For given scalar field potentials with n = 1/2, 2/3,1, 2, 3,4, 5, 6, A/^ = 60 and ( = 10 -8 , the non-minimally kinetic 
coupled DGP model lies inside the 95% CL Planck+WP+highL data for n = 1/2, 2/3,1,2, and lies inside the 
95% CL Planck+WP+highL+BICEP2 data, for n = 1,2, 3,4, 5, 6 (see Fig.9 and table III). 

• For given scalar field potentials with n = 1/2, 2/3,1, 2, 3,4, 5, 6, N = 60 and ( = 10 -10 , the non-minimally 
kinetic coupled DGP model lies inside the 95% CL Planck+WP+highL data for n = 1/2, 2/3,1,2, and lies 
inside the 95% CL Planck+WP+highL+BICEP2 data, for n = 1,2, 3,4, 5,6 (see Fig.10 and table IV). 

In conclusion, in the study of inflation using the non-minimally kinetic coupled DGP model, we found that among 
the suggested potentials and coupling constants, subject to the e-folding N = 60 required by inflationary scenario, 
the potentials V (</>) = cr^, V (</>) = crcj) 2 and V (</>) = acj) 3 provide the best fits with both Planck-fWP+highL data and 
Planck-f-WP+highL+BICEP2 data. 
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